N94-34197 


parameter identification of a rotor supported in a pressurized 

BEARING LUBRICATED WITH WATER 

J° hn W. Grant, Agnes Muszynska, and Donald E. Bently 
Bently Rotor Dynamics Research Corporation 
Minden, Nevada 


Y 


abstract 

Applying a nfntynltono^ Water ~' U ^ ricatcd bea ""g was developed 

amphtude as an input, rotor vibration respond dat* 8 perturbation with a con.tan 

format,. Using this data, the parameters ” the ZT “ B ° de and Dynamic Stiffnes 

rotor/bearing model was represented by the «n 1 !° t ° r / bea "“8 Wtem were identified. Th 

mode ’ wi,h *• <*««« or l nuw SuSizr 1 "* of ,h ' f, ' , ‘ ut "" 

1- INTRODUCTION 

structures. ^pplicS of Lda^leibg'tt rtmcturf ^ Pr ° bl ' mS ° f mechan j ca 

Vstem* parlme^e^^lch^ can in ^ the^vduat^^t” data '^.'^n^rionTih! 

m numerical modelling or engineer assist (expert svstemT ^ l ° ^ 8tablbty of the 8 y«tem, and 
of perturbation testing and parameter identification Pr ° gram “' Thi " pa P« Presents results 
externally pressurised water-lubricated bearing Th * 8Upported in one rigid and one 
system identification was developed in references £^ 11 ? 't™* 1 *'* r " odel «“d for the rotor 
mput force to measured output response (henceforth lrn ^ !" etbod of relating the applied 
«rve. very well for the “ the dynamic stiffne8s technique) 

experimental testing of the rotor/bearing sy,!m This *T ***■ d ' Vel ° ped f ° r the 
amphtude perturbator used to apply the force m thl' ^ °? C ° n81st8 of a constant force 

±3000 rpm (perturbation forward and backward! The" t<!8t ° VCr a fre< 3 uenc y range of 

the force has a constant amplitude, independent af fr ^ ° f th ‘ 8 perturbator “ that 

Perturbation. HI*. «„. b l«. to'.w”. mv ^ •"» 

the low frequency range, therefore it offers mn » ^ (high signal to noise ratio) data in 

tutor, with fluid intention. ’ ™ m °" ■<‘»M«ti„n of th, low,., mod,, of 


2 . 


experimental test rig 


The experimental rotor is shown in Figure 1 The . f 

carrying a balanced, concentrated mass^t the brasa-sf* ^ ^f 818 ^ of a 8tabIesa «teel sha 
diameter and 10 mil diametrical clearance inside th iT ? Ved ’ aIum,num journal of 1.990 inc 
by a relatively rigid, bronze (oilite) bearfog and hi otte t' ** “ 8UPP ° rted on «•« - 
bearing This cylindrical bearing has fouf’equall! BD ‘ ed * hydro ‘ ta * 

Figure 2. These pocket, are fed through rXlO iTs.n’h^ 0 ^^ P ° ckeU ' “ I 
contains, also, separately controlled tangential , *• ' h dlameter capillaries. The bearin 

,h ' - ™ 1 - «* mJ r£ [ n 4. b ' h0 "'™ 
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Fie 1 Experimental rotor rig with constant force perturbator, A) 0.5 hp electric motor', B) P* 
Flg - 1 conUol probc, C) Keyph-r® prob,, D) flexible coupling, E) bronsc sliding b-mg,J) 
shaft, G) slim disk for balancing, H) spring support (with stiffness 38 lb/in), ) ) ? 

current displacement probes, J) experimental pressumed bearing, K) constant force 

perturbator, L) optical Keyphasor. 

antiswirl port 
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Fig. 2 Externally pressurised bearing with antiswirl ports. 


Utera? 1 vibration 'data, on, pair of XY ,dd, ,o„,nt di*lac,m,». prob« wa, nroonUd to 
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was uwd to ttlr t0 thC j ° Urnal ' A com P uterizcd acquUition and processing system 

syster/consists of a rt^stee^A^coMecST tT^ *° Pe . rtU . r ^. the rotor /bcaring system. The 
a sealed rolltog element bearing The shaft ° * JOUrn f^ Wltbln l ^ e hydrostatic bearing, using 

b«m g . A cLunt rlXli™ .X’h«Xt °Y* ° lh " '"i* » (oi»4 

rolling element bearing (Fig 31 The .atari . C next to l* 16 journal through a 

hp bidirectional electa molo/ I 'utt Y""? ?' ‘ h ' pCTlurb * ,i ‘ > “ »“ drived from a 0.5 
and acceleration. The mto retlnel 7*5. '°‘" r ™ to “" ,rol P«»urUtion .peed 

perturbation frequency 'JwZXX ™ '° P, “ T.‘ Urb f ' fot “' »“ Altered to the 
•- « b — ed a reflective SrrpSX *° “ o^* 1 



I 

Fig. 3 Constant force perturbator (not to scale; eccentricity £ exaggerated). 

3-, CONSTANT FORCE AMPLITUDE PERTURBATOR 

?n~bltZl^f tb *‘“" f ” re ' «• ‘ b < — wa, generated 

amplitude becomes very small for low Dertu^bat^'f 18 prop0rtlon A aI to rotatlvc speed squared, its 
Provides a constant force t ^ T de ‘ igncd that 

range. The principle of the ooeration JM? T the data accurac y m the low frequency 
excitation, known in vibration thTorv- a calib^l ?"! perturbator “ simila r to the kinematic 
across a known radial stiffness (K) 7{ the ^ ] ! ng . * c . ctor of displacement (£) is exerted 

known rotating force vXThat L ann^ within the P«t^bator. This create, a 

obtained by adjusting the eccentricitv of an°' * rot °f und * r test - Tbe known displacement is 
obtain the rotating force 8Wee n f " Cahbrated *> «hown in Figure 3. To 

acceleration to a direction the same ^Tth^' * p * rturbatl ° n disk is rotated with very low 

the direction (reverse perturbation) The eto^mat' ^ perturbation ). 0f opposite 

Corning® 3112 RTV silicon rubber The .*■«- ™ aten ^ U8ed for the perturbator was Dow 
separately identified by the use of’ dvnaiJr^fT’ d " nPI . ng ' inertia of this ma terial were 
found that maximum inertia ***** ^ 14 *« 

* -1 compared .o .he luiee, .heySe %£ tt&tSZXT * 
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4. WATER DELIVERY SYSTEM 

To study hydrostatic bearing/seals, a high pressure water system (1000 psi O 10 gpm) was 
developed. It consists of a 50 gallon water storage drum, low-pressure feed pump 100 psi), 
high-pressure multistage centrifugal pump (400 psi) driven by a 5 hp 3 phase 230 volt electnc 
motor rotating at 3450 rpm. This delivers water through a pipeline which contains 80 me , 
micron, and 1 micron filters, and several gauges and valves. Water is returned to the 50 gallon 
storage drum using a sump pump. Water temperature is openly controlled y using a 
submersion heater in the storage tank, with cooling provided by a fan blowing air at the storage 

tank surface. 

The water is delivered to the hydrostatic bearing using three specially designed distribution 
manifolds. One manifold leads to the four radial ports, and the other two to the eight anti.wirl 
ports of the test bearing. Each port is fed from the manifold to the port through » manuaHy 
controlled flow valve, and a high-pressure hose. Each of the high-pressure feed lines « 
equipped with pressure gauges to measure the pressure drop at the bearing port. Pressure is also 

monitored at the manifold. 

Flow is measured using the sump pump. The sump holds a calibrated 2.5 gallons before 
pumping, thus by timing this volume, the total flow rate is determined. Such calculation is on y 
valid when testing one set of ports at a time, and when the journal rotates concentrically within 

the bearing. 


5. MATHEMATICAL MODEL 

The mathematical model of the fluid dynamic force and the corresponding dynamic stiffness 
components of the rotor/bearing system were developed in [1,3]. 

The model of the considered system is as follows: 

Ms + Mf(s - 2jAQs — A 2 H 2 *) + D(s— jAfls) + (K 8 +Kb+K)s = K£e^ ^ , j=\/-i 
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m: “-j - rt k * “ d m « 

radial stiffness and damping respectively Jd i ^ ' K ® ^ D the fluid 

ratio (the rate at which «,f ^ * ‘ 8 th * flmd average circumferential velocity 

t ne rate at which the fluid damping force rotates), with u and Q beimr Derturh*tmn 

tt« P*,to. grated eccentricity garnet*., itl v „« „ d ^ 

“;“!”.'“>“"“.“« = Ae J <“ + “).he re A and a are amplitude and pk« 
■eepecttvely. The follomng rel.l,on.hip. «re further u«d in the identification procedure, 


Direct Dynamic Stiffness = K g +K B +K— Mf(w— Afl) a - Mw a = M CO s (6-a) 
Quadrature Dynamic Stiffness = D(a>-An) = — sin(5-a) 


( 1 ) 
( 2 ) 

(’l) y «tnd m ('J) to'det™t“t® ZZ fTJh h ^,™Z* l ' d 10 th ' “ ** 

6. EXPERIMENTAL TEST RESULTS 

!° de ! crmlne the cffecU of pressurised water input through the 
wTlpp IZlt fv 7^“^ ( Tr •? P ° rt8 8hut off >’ The P" 88 ™ at the radial port mfnifcdd 

constiSt at 20 p f The w^er !] t u° P T*? thr0Ugh the was maint ^ned 

consecutive test? n • f fl ° W thr ? ugh the bcarin g v “ied from 3.9 to 4.6 gal/min for four 

and 3000 rpm respectively. ^ C ° n8CCUtlVe tcst8 ’ the rotativc 8 P ccd wa8 constant 0, 1000, 2000, 

£K = 1 535 lbl'at^ — SlO^Th* ^ & rotatmg 8Wce P frequency perturbation force vector of 
ru J A ” 310 ’ Thc rcs P° n8€8 to perturbation forward (+u) and backward (-u) 

tTce tJatTh he perturb . atic ; n ftequende. and are presented in the B^efo" a t(Fi^ ( No^ 
«l t i theTr ampUtude V“ thc ~™ d - d ™ge of frequencies are rather flaVnot the 

EL fa’ ' a ; t i nbalanCC f0ICC - N ° te - ak °' that thcre “ d ^ 

r«rS “ “d Phases around the area when the perturbation and rotation frequency 
system. ue to an interference of the rotative frequency vibrations in the filtering 

Figure 6 presents the dynamic stiffness plots that were obtained using Eqs. (1) and (2) and the 
m or response shown in Figure 5. Each plot consists of 570 data points ! To identify tt 

cLffidentTtf a hi?h rametCra ’ S' Pl ? 8 W " C bC8t flttcd by P ol y nom *als. Interesting enough, the 
Stiff* b gber P° wer than the *econd for the direct stiffness, and first for the quadrature 

“ XSTJS T f ' 8mallerthaD th °“ of powers. Therefore' a ptatoU 

r ,h ‘ pku ° r pieu " e - Tt ' «■— ** « 

,h< “ of lhe dynamic stifTnese model (1), (2), and cou- 
m Table 2. ? “ g “ “ fre 1 u ‘” c y unlt *i ‘he eyelem paiamelete can be identified, aa ehoven 

^The^tLot' r^t'^d jrT 1 !?.? r To'obto“ !“ ““ u t Ct 

for eero epeed, the modal maaa wae obtained by averaging the valuee at all other rotatile a'peede* 
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Pig. 5 Rotor vertical response (filtered to perturbation frequency) to a forward and backward 
perturbation force rotating vector of K£ = 1-535 lb. at 6 = -310 , with the 
water-lubricated bearing radial port pocket pressure of 20 psi and a rotative speed of (a) 0 
rpm, (b) 1000 rpm, (c) 2000 rpm, and (d) 3000 rpm. The horisontal responses which 
were very similar are not shown. Response amplitude measured peak to peak. 
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Table 1. Best Fit Coefficients of the Polynomials in Fig. 6 


Rotative 
Speed Q 

Direct Dynamic Stiffness 
A + Bm + Cm 2 

Quadrature Dynamic 
Stiffness Em + F 

rpm 

A 

flxl0~ 3 

CxlO -5 

E 

F 

0 

1182 

2.126 

-8.084 

0.400 

3.424 

1000 

1226 

49.30 

-8.557 

0.420 

-150.5 

2000 

1205 

113.6 

-8.363 

0.408 

-318.9 

3000 


152.0 

-7.707 

0.402 

-530.2 


Table 2. Identified Parameters of the Rotor/Bearing System 


Rotative 
Speed, Q 

D 

n 

M f 

K+Ks + Kb 

M 


rad 


kg 


lb x sec 2 

kg 

lb 

N 

lb x sec 2 

kg 

sec 

m 

sec 


inch 

xlO -3 

inch 

m 

xlO 3 

X 

S . 1 

O 

0 

■ 

669 

— 

5.94 

1.04 

1182 

207 

1.43* 

0.251* 

1000 1 

105 

4.01 

702 

0.358 

6.28 

1.10 

1235 

216 

1.52 

0.267 

2000 

209 

IB129 

E31 


6.64 

1.16 

1249 

219 

0.99 

0.173 

3000 

314 

■ana 

WSSM 

IEES3 

5.25 

0.92 

mmm 

231 

1.78 

0.312 


♦Average of the column 


The average mass, M a vg = 1.43 * 10"* lb sec*/inch (0.251 kg), was then subtracted from the 
second order coefficient C , providing the fluid inertia effect value at 0 rpm. 

As it was mentioned in reference [2], there is a possibility that the fluid radial damping and fluid 
inertia forces rotate at different rates of A . A method to determine the fluid inertia rotation 
rate, Af, is to solve the equation resulting from equating Eq. (1) with the best Fit coefficients and 
solve it for MfAf (see Table 1 and Eq. (1)): 



where B is the first order coefficient of the direct stiffness best fit parabolas (Table 1). If the 
average model mass, M aV g» i» used again, the fluid inertia effect can be calculated from the 
second order coefficient of the parabola: 



where C is the second order coefficient of the direct dynamic stiffness best fit parabola (Table 
1). The ratio Af is calculated by dividing Eq. (3) by Eq. (4). The results are summarised in 
Table 3. They show that the values A and Af are, actually, very close (maximum difference 
is 7.2%). 
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Table 3. Identification of the Fluid Inertia Average Circumferential Velocity Ratio 



Figure 7 present* the results from another set of perturbation tests in which the rotative speed 
varied, while the perturbation force was constant with Ke = 1.535 lbs and 6 = —310* . The 
quadrature dynamic stiffness data were identified as parabolas or straight lines by best fitting 
polynomials (Table 4). The corresponding quadrature stiffness components D and A that 
were calculated from the data are summarized in Table 5. Generally, D resulting from these 
cs is s g y lgher, and A slightly lower than obtained from the previous tests. 

All test results were not compensated for the water temperature changes, which, within each test, 
varied from 0.3 F to 0.8 F within the range of 90.4* F to 96" F. The data acquisition system 
was set on the average value of the proximity transducer sensitivity, 0.2 v/mil. Due to slightly 
increasing temperature during the tests, the transducer output sensitivity varied from 0.198 v/mil 
to 0.199 v/md thus the amplitude readings were 0.5% to 1% higher than the actual ones, and 
consequently, the dynamic stiffness values became from 0.5% to 1% lower 
1 500 -t 


, , ». » » cj= 500 rp[m 

j j 00000 10 >= 1000 rpm 

l __ L500 _rpm_ 

*>*>1 j io>=2005 rpm 

I kj=2500 rpm 
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V I KJKJKJ 5000 
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Fig. 7 Quadrature dynamic stiffness versus rotative speed for constant values of perturbation 
trequency. The solid lines are best fits using second order polynomials. 
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Tabic 4- Best Fit Coefficients for the Polynomials of the Quadrature Dynamic Stiffness From 
Fig. 6 Versus Rotative Speed. 


Perturbation 
Speed ©, rpm 

Second Order Fit 
G + HCl+I n 2 

First Order Fit 

JQ + L 


G 

tfxlO -3 

JxlO -6 

JxlO -3 

Pi 

500 

216 

-136 

-13 

-207 

287 

1000 

443 

-123 

-19 

-228 

543 

1500 

573 

-115 1 

-17 

OO 

o 

1 

662 

2000 

835 

-159 

-13 , 

-231 

904 

2500 

1097 

-188 

-8 

-229 

1137 


Table 5. Identified Parameters for Quadrature Dynamic Stiffness Using Eq. (2) and 
Polynomial Fits From Table 4. 


Perturbation 

Second Order Fit 



1 

Speed, © 

Parameters 

First Order Fit Parameters 1 


D 

H 


■ 

rpm 

rad 

lb * sec 

kg 


lb* sec 



sec 

inch 

sec 


inch 

sec 


500 

52 

4.13 

723 

0.3 15 

5.47 1 

958 

0.361 

imm 

105 

4.23 

741 

0.278 

518 

907 

0.420 

1500 

157 

3.64 

637 

0.301 

4.22 

H£l 

0.471 

2000 

209 

3.99 

699 

0.381 

4.31 


WEsm 

2500 

262 

4.19 

734 

0.429 

4.34 

760 

| 0.504 | 


The fluid radial damping calculated using the well-known formula: D = 7CT}r(l/c)* (where 7} is 
fluid dynamic viscosity, r, 4 and c are bearing radius, length, and radial clearance respective- 
ly) provides the values from 3.45 lb sec/inch to 3.22 lb sec/inch for water temperatures from 
90.4* F to 96.5* F and corresponding viscosities varying from 0.7577 c poise to 0.7076 c poise 
through the entire set of tests. The identified values of D were, therefore, lower than the calcu- 
lated ones. The proper temperature compensation of the results would decrease the difference. 

7. CONCLUSIONS 

Developed over two years ago, the high-pressure water rotor rig can be used for dynamic testing 
of bearing and seals. This paper presented test results of an externally pressurised bearing for 
only one pressure value. Applying the nonsynchronous perturbation testing which has been 
extensively used by the authors through the last 12 years, the fluid force parameters were 
identified. The results confirm, again, the meaningfulness of the fluid circumferential average 
velocity ratio as the important element of the fluid force model in lightly loaded lubricated 
bearings and seals. The results also indicate the existence of the fluid inertia effect which is 
relatively high, about five times larger than the rotor modal mass. 
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The use of the constant force perturbator has enabled the acquisition of very clean vibration data 
(higher signal, thus better signal to noise ratio) in the low frequency range, and hence clean 
dynamic stiffness data was produced. This data makes the identification procedure more 
accurate, and allows for advancements in adjustments of the mathematical model of the 
rotor/bearing system. 


NOMENCLATURE 


A t B t C t E t F,G t H t I,J t L 
A, a 

D^Mf 

K S ,M 

K 

*=x+jy 

6,e 

A 

Af 

OJ 

n 


Best fit polynomial coefficients 

Rotor response amplitude and phase respectively 

Fluid radial damping, radial stiffness, and inertia effect respectively 


Rotor modal stiffness and mass respectively 
Perturbator radial stiffness 

Rotor lateral (x=horizontal, y=vertical) displacement 

Perturbator angular orientation and eccentricity respectively 

Fluid circumferential average velocity ratio at which fluid damping force 

rotates 

Ratio of the fluid inertia force rotation 
Perturbation frequency 
Rotative speed 
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